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Abstract. In NMR-based quantum computing, it is known that the controlled- 
NOT gate can be implemented by applying a low-power, monochromatic radio- 
frequency field to one peak of a doublet in a weakly-coupled two-spin system. 
This is known in NMR spectroscopy as Pound-Overhauser double resonance. The 
"transition" Hamiltonian that has been associated with this procedure is however 
only an approximation, which ignores off-resonance effects and does not correctly 
predict the associated phase factors. In this paper, the exact effective Hamiltonian 
for evolution of the spins' state in a rotating frame is derived, both under irradiation 
of a single peak (on-transition) as well as between the peaks of the doublet (on- 
resonance) . The accuracy of these effective Hamiltonians is validated by comparing the 
observable product operator components of the density matrix obtained by simulation 
to those obtained by fitting the corresponding experiments. It is further shown that 
an on-resonance field yields a new implementation of the controlled-NOT gate up to 
phase factors, wherein the field converts the 1^ state into the antiphase state 2I^lf, 
which is then converted into the desired two-spin order 21^1^ by a broadband n/2 
pulse selective for the A spin. In the on-transition case, it is explained that while 
a controlled-NOT gate is approximately obtained whenever the radio-frequency field 
power is low compared to the spin-spin coupling, at certain specific power levels an 
exact implementation is obtained up to phase factors. For both these implementations, 
the phase factors are derived exactly, enabling them to be corrected. In Appendices, 
the on-resonance Hamiltonian is analytically diagonalized, and proofs are given that, in 
the weak-coupling approximation, off-resonance effects can be neglected whenever the 
radio-frequency field power is small compared to the difference in resonance frequencies 
of the two spins. 
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1. INTRODUCTION 



The c-NOT (controlled-NOT) gate is of central importance in quantum computing 
( [Barenco et al.| |1995a| , |Barenco et al.| |1995b ). It is well-known ( [Lloyd 1993 ) that in 
quantum computers based on frequency addressing of their q-bits, such as solution-state 
NMR spectroscopy on spin | nuclei (|Cory et al.||1997| , |Gershenfeld and Chuang||T997|) , the 
c-NOT gate can be implemented by the application of a two q-bit transition Hamiltonian 
of the formFj] 
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where E± = |(1±2J B ) is idempotent (cf. ( |Hatanaka and Yannoni| [1981|) ). Using NMR 



spectroscopy ( |Cory et al.| |1998|) , we have previously demonstrated that irradiating 
exactly one peak of a doublet in a two-spin system with a power uj\ much less than 
the coupling between the spins 27rJ AB transforms the equilibrium state in accord with 
the corresponding propagator, i.e. exp(— LtH trn ) 
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In NMR spectroscopy, this is often called Pound-Overhauser double resonance flSlichtci 
1990|) . The nonzero phase factors in this matrix can be equalized by a 7r/2 evolution 
exp(— l(-k/2)I b ). Nevertheless, the application of this procedure to a (mixture of) 
superposition states quickly shows that the effective Hamiltonian is not so simple. The 
transition Hamiltonian is only a nonphysical approximation, which cannot predict all 
the details of the spins' evolution under selective irradiation. 

In this paper the effective Hamiltonian in a rotating frame ( [Ernst et al.| |1987 . 



Slichterj |1990|) is derived, which fully describes the evolution of a weakly-coupled, two- 
spin system under a monochromatic RF (radio-frequency) field, both on-resonance as 
well as on a single transition. A pictorial representation in terms of effective fields is 
described, which provides an intuitive description of the spin dynamics under these 
Hamiltonians. The results of NMR experiments are presented, which demonstrate 
that the superposition states that evolve under monochromatic RF fields are consistent 
with those obtained from simulations using these effective Hamiltonians, and with the 

B | and that off- resonance 



effective fields picture. Assuming u\ <C 2n\J \ <C \uj 
effects can be neglected, it is proven that the transition and effective on-transition 

f Throughout this paper we shall be making use of the product operator notation that is widely used 
in NMR spectroscopy (Sorensen et al. 1983), and the rules for manipulating these symbols derived from 
geometric algebra ( Somaroo et al. 1998| ). 
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Hamiltonians are equivalent, in the sense that the corresponding propagators are 
approximately equal up to conditional phases (i.e. a diagonal matrix of phase factors). 
This implementation of the c-NOT is shown to be exact (up to conditional phases) 
for u\ = 27rJ AB /\/4n 2 — 1 where n > is an integer. It is further proven that the 
application of an on-resonance field with uii = ti~\J ab \ followed by a broadband n/2 
"soft pulse" (covering the entire doublet), both on the A spin, likewise implements the 
c-NOT gate up to conditional phases. In both cases, the conditional phases are derived 
explicitly. In appendices, the on-resonance Hamiltonian is analytically diagonalized, 
and the assumption that off-resonance effects are negligible whenever uoi <C \uo A — uo B \ 
is rigorously justified under the weak-coupling approximation. 



2. THE POUND-OVERHAUSER EFFECTIVE HAMILTONIAN 

2.1. Derivation of the effective Hamiltonians 

Given a Hamiltonian H = H(t) and unitary transformation U = exp(— iGt) (where 
G = G is Hermitian), the evolution of the transformed density matrix p' = UpU is 
given by QErnst et al.| [T987| , ISlichterj [1990D 



p = UpU + UpU + UpU 

= - iGp' + i[p', H'] + ip'G (3) 
= l[p',H' + G] = i[p',H eS ] . 

where H' = UHU and -ff e fr = H' + G. The weak-coupling Hamiltonian of a two-spin 
system is| 

riAB , ,AtA , RrB , tAB rA tB t a\ 

H = - uj 1 z - u Q l z + 2ti J l z l z , (4) 

where u A , u B are the Larmour precession frequencies of the spins, and J AB the scalar 
coupling between them. The applied RF field Hamiltonian has the form 

H RF = u ie - tG (I A + I*)e- tG , (5) 

where the RF power transmitted to the spins is (assuming for convenience that the 
system is homonuclear) iO\ > 0, G = ^{I^ + ^z)j an d ^2 is the frequency of the RF 
field. 

These results show that the time-dependence can be removed from the Hamiltonian 
H = H AB + i?RF by transformation to a frame rotating at the frequency ui2. If the 
frequency u 2 = uj a + nJ AB matches the | 01 ) <-> 1 11 ) component of the A-spin doublet, 
this yields the time-independent "effective Hamiltonian" 

H cS = ttJ ab (/ z a + 2/ A / z B ) + (u; A - uj B + ttJ ab )/ b + Wl (J x A + J B ) . (6) 

X In quantum computing, the ground state is generally indexed by "0" , whereas NMR spectroscopists 
typically put the spin "up" (parallel the field) state before the "down" in their matrices. These two 
conventions agree only when the gyromagnetic ratio is positive, as will be assumed in this paper. 
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Figure 1. Effective fields picture of the evolution of the irradiated A-spins in the 
on-resonance case. There are two subpopulations of molecules: that in which the B- 
spin is "up", and that in which the B-spin is "down". In these two subpopulations, 
the magnetization due to the A-spins is initially aligned with the applied field (broad 
banded up arrow). The RF field power is lo\ = tt\ J ab |, which in the co-rotating frame 
yields the effective fields shown with the thin solid arrows inclined at 7r/4 from the 
z-axis. After a time t = tt/((uji) 2 + (tt J ab ) 2 ) -1 / 2 , these fields have rotated the A-spin 
in both subpopulations by an angle of 7r to the ±£-axis, which corresponds to the 
antiphase state I^I^- 

In the event that the RF field is placed on-resonance, u>2 = uJq, and hence 

H cS = 2vrJ AB ^/ z B + K A - u*)I B z + + J B ) . (7) 

Note that since G = ^2(1? + / B ) commutes with / A / B + I y Iy, if we use the full 
strong coupling Hamiltonian 

ff AB = - ^ A / A - Uq -T b + 2vr J AB (I X A J B + J A J y B + J A J Z B ) . (8) 

instead of its weak coupling approximation (Q), Eqs. @ and (^) need only be modified 
by the addition of the term 2ttJ ab (J a J b + I A I B ). 
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Figure 2. Effective fieids picture of the evolution of the irradiated A-spins in the 
on-transition case. As in Figure [l], the magnetization due to the A-spins in the two 
subpopulations of molecules is initially aligned with the applied field (broad banded 
up arrow). The RF field power is wi <C 7r| J AB |, which in the co-rotating frame yields 
the effective fields shown with the thin solid arrows. The one corresponding to the 
irradiated transition is along the x-axis, while the other one is displaced towards 
the z-axis by an amount 2irJ AB . After a time t = tt/ui, the A-spin in the first 
subpopulation has been rotated by tt to the — z-axis. If u>\ is set so that at this time 
the A-spin in the other subpopulation has rotated by a multiple of 27r, this yields an 
exact implementation of the c-NOT up to conditional phases. 



2.2. The effective fields picture of the evolution 

A physical picture of the spin dynamics may be obtained by breaking the NMR spectrum 
up into its individual resonance lines, and specifying a distinct effective field for each. 
This picture assumes we are dealing with an equilibrium density matrix in the high- 
temperature approximation, so that there are no spin-spin correlations, and neglects 
relaxation effects. In a frame rotating at the transmitter frequency ui2, the residual 
static magnetic field at each resonance ujq is equal to the frequency offset ujq — uj 2 . Each 
of these residual fields may be identified with a subpopulation of the molecules present, 
wherein the other spins are aligned so that their couplings with the spin in question 
cause it to resonate at the frequency of the residual field. The applied radio-frequency 
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field contributes a transverse component of strength iO\ to the net effective field for each 
resonance. With care, all the results of this paper can be derived by purely geometric 
means from this "effective fields picture" . 

In the on-resonance case with uj\ = tt\ J ab |, the effective fields for each resonance 
line of the A-spin doublet form angles of ±7r/4 with the z-axis, as shown in Figure 
III. Thus the equilibrium magnetization components for each resonance counter-rotate 
from their initial positions along the z-axis to opposite directions along the x-axis after 
a time "t" given by it = t\J (ooi) 2 + (it J AB ) 2 . This vector configuration represents the 
antiphase state I A I^, which may be converted to the I A I^ state expected after a c- 
NOT gate by a broadband 7r/2 |/-pulse selective for the A-spin. In the on-transition case, 
in contrast, the on-resonance component of the magnetization experiences an effective 
field in the transverse plane, while the other component experiences a field which is, 
for uj\ <C \u A — (J® I, very nearly along the z-axis. Since it is exactly along the z-axis 
only in the limit uoi — > 0, however, for any uj\ > it nutates away from the z-axis, as 
shown in Figure |2[ An exact implementation of the c-NOT gate is nevertheless obtained 
when Ui is chosen so that this component is rotated by an integer multiple of 2n back 
to the z-axis in the time it takes the other component to rotate to the — z-axis, namely 
t = 

3. EXPERIMENTAL VALIDATION OF THE EFFECTIVE 
HAMILTONIAN 

In order to demonstrate that these effective Hamiltonians correctly describe the 
evolution of a weakly-coupled two-spin system under a monochromatic RF field, we 
used a solution of triply-labeled 13 C alanine (C0 2 — CH[— NH^] — CH3) in D2O. The 
carboxyl carbon was treated as the "A" (target) spin, the alpha-carbon as the "B" 
(control) spin. Since the carboxyl to methyl carbon coupling constant of ca. 1.4 Hz. is 
less than the peak width, its effects could be ignored, and the effects of the protons 
eliminated by decoupling. The following three experiments were performed at 16 evenly 
spaced time points in the interval from t max /16 to t max : 

(i) On the u A resonance with uj x = ir\ J AB | (t max = \/2/\J kB \). 
(ii) On the uj a + vrJ AB transition with u x = 2vr| J AB | (t max = 1/| J AB |). 
(iii) On the u A + nJ AB transition with uji = 4n = 2n x 2 Hz. (t max = tt/o^x = § sec). 

Experiments (i) and (iii) validate the correctness of the effective Hamiltonian in the 
two cases that may be used to implement a c-NOT gate, while (ii) validates the theory 
for case in which the evolution of coherence is more complicated. All the experiments 
were performed with the spin system initially at equilibrium in a 9.4 Tesla field. 

After the RF field had been applied for each of the 16 time periods used, a 24k 
point FID (free- induction decay) was collected over 590 msec, zero-filled to 48k, and 
Fourier transformed to yield the complete spectrum. A 128 point window centered on 
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Figure 3. Plots of observable product operator components of the density matrix 
versus 16 equally spaced time points (in percent of i max ) for the on-resonance case 
(i) with uji = 7r|J AB | and t max = V2/\J AB \. The symbols in this plot (and the 
corresponding product operator components) are "o" (Ix)> " D " Cy)i "A" (-^x-^z)i 
and "O" (/ A if). The solid lines are the simulated results for these components, while 
the dashed line is the simulated evolution of the diagonal I A component. 



the peaks from the A-spin was taken and the absorptive in-phase (I*), absorptive anti- 
phase (21^ if), dispersive anti-phase (21^ if) and dispersive in-phase (Jy ) components 
of the density matrix extracted from it. This was done by a least-squares fit of a linear 
combination of model peak shapes for each of these four components to the spectrum 
window, where the model peak shapes were computed using a coupling constant of 54 
Hz. together with a Lorentzian peak shape of half-width 0.85 Hz. 

In order to compare these results with the theory, the results predicted from the 
above effective Hamiltonians were also computed by numerical simulation, using the 
same coupling constant and resonance frequencies. These results were normalized so 
that the initial 1^ state had unit norm, and the experimental fits scaled so that the root- 
mean-square values of the experimental results were the same as the simulated results. 
If necessary, the spectrum was adjusted by means of a first-order phase correction so as 
to obtain the best possible fit, as judged visually. The simulations and final fits to the 
data for each of these three series of experiments are shown in Figures ^| - 15|. It should 
be clearly understood that the simulated curves were not fitted to the experiments save 
by scaling. 
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Figure 4. Plots of observable product operator components of the density matrix 
versus 16 equally spaced time points (in percent of i max ) for the on-transition case 
(ii) with lui — 2n\ J AB | and i max = 1/| J AB |. The solid lines are the simulated results 
for these components, while the dashed line depicts the evolution of the diagonal J A 
(starting from 1) and / A / B (starting from 0) components. The symbols in the plot are 
the same as in Figure ||. 



Given the many sources of systematic error present in NMR spectroscopy ( [Hoch aiicl 
Sternj |1996| ), and the fact that the fits were not systematically optimized with respect 



to the nonlinear parameters (i.e. the coupling constant, resonance frequencies, peak 
widths/shapes and spectrum phases), the match of the experimental to the simulated 
results is strong evidence for the validity of the theory. The corresponding simulations 
using the transition Hamiltonian (not shown) are essentially the same, except that the 
high-frequency oscillations are not present. 



4. DERIVATION OF THE CONDITIONAL PHASE FACTORS 

4-1. The on-resonance case 

Under the assumption of weak coupling, it is shown in Appendix B that if uj\ <C 
|u; A — uj b \, the off-resonance effects due to the I B term in H c k can be ignored. This 
can be understood intuitively through the effective fields picture for the B, since under 
these fields are essentially along the z-axis and hence only induce phase shifts. Thus let 




Figure 5. Plots of observable product operator components of the density matrix 
versus 16 equally spaced time points (in percent of i max ) for the on-transition case 
(iii) with uj\ — An and i max = \ sec. The solid lines are the simulated results for these 
components, while the dashed line depicts the evolution of the diagonal I z (starting 
from 1) and (starting from 0) components. The symbols in the plot are the same 

as in Figure ||. 



ci = 7rJ AB /u;i and c 2 = (uJq - cj b )/u x , (9) 



so that 

HVo* = c 2 / z B + (I A + 2 Cl I*I*)(E B + El 



, 2 f B + (J A - Cl J A )£ B + (JA + Cl /A)«B ( U) ' 



Since all three terms in this expression commute, it follows that exp(— UH cS ) 
= e -«.o; 1 t(zA- cl 4)EB e _ 1 ^ lt( jA +clJ A )E B e _ la;ite2l B 



where by the formula for the exponential of an operator multiplied by a commuting 



idempotent ( Somaroo etaLl 11991) , exp(-iwi* (J A ± cJ^El 



E% + fcos(iu;it > /i+^ 



± Cl J A ) 



sin(icJit v /TT^) J -E 



(12) 
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For u)i = n\J AB \ and t = n / (y/2ui) , this implies exp(— LtH® s ) 

= - iV2 ((I A - I k )E B + (I A + I A )£ B ) e-^l^ 2 

= - 2l (l A E B + I A E B + ) e ™$/* e -™»$/^ (13) 

_ g-Mrl£/2 e -cnI*E B _ e -wI*E\ e -LTTC 2 I*/V2 



where have we used the fact that J y anticommutes with I A and J z to move its 
exponential to the far left. Thus the applications of an on-resonance pulse with 
uj\ = 7r|J AB | and t = n/(y2ui), followed by an ordinary soft — %/2 y-pulse on spin A, 
yields the c-NOT gate up to a conditional phase prefactor of exp(— t7r(/ A /2 + (c 2 v / 2 — 
l)J z /2 + I A I B )). By moving exp(t7r/ A /2) to the far right instead, one sees that one 
could just as well precede the on-resonance pulse with a soft +7r/2 y-pulse on spin A. 



4-2. The on-transition case 

We begin again by breaking H° eE = H eS — u>il x into three commuting parts as above, 
namely 



(7rM)i*° ff = n((c 1 + c 2 )I B + 2c 1 I A E B + + I A (E B + + E B )) 

2 )/ z B + (2ciJ 
P + Q + E B +Q E B 



n (( Cl + c 2 )/ z B + (2 Cl J A + I A )E B + + I A E B ) (14) 



where 



P = tt( Ci + c 2 )/ z B , Q + = 7t(2 Cl I A + 1$), and Q = nI A . (15) 



Note that P is diagonal (i.e. along the z-axis) and that H tIIi = Q E B , but since Q + E B _ 
is not diagonal, exp(— b{i: / uj\)H eS ) and exp(— iH tvn ) do not simply differ by conditional 
phases as in the on-resonance case. Nevertheless, Q + E B _ can be readily diagonalized, 
because 

Q + = 7r( e -^/ A e ^ A )v^T4^, (16) 

where 9 is given by arctan(l/(2ci)). 

To show that the effect of this diagonalization upon the propagator exp(— ttH° cS ) 
is small, we evaluate the norm of the difference of the propagator with and without 
this transformation.^ Since I A and Q = ttI a don't commute, we operate with 

§ The norm we use is the square-root of the geometric algebra "scalar part" (denoted by angular 
brackets (•)) (| omaroo et al.| n99q) of the product of the quantity with its Hermitian conjugate. For 



the two spin system considered here the scalar part is four times the trace in the usual product Pauli 
matrix representation, and hence our norm is just twice the standard Frobenius norm. 
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exp(— iQIyE^) instead of exp(— iOI A ), obtaining 



e 

2 



= 2-2 (^ e -'-^+ I ^ EB + e ^V^4^ E i 

// 2 Cl / A + / A \ ( 17 ) 

= 1 ~ (^cos(f - i i 1 ^_ ^sm(|yi^)j 

cos(f v^+^f) + t2I A sin(f V /I^) 

= l-cos 2 (f v ^f)-sin 2 (f v ^)^= 

V 1+4c i 

= sin 2 (f V /I^) (1 -2c 1 / v ^f) 
The second factor can be converted into the simple bound 

l + 4c 2 -2lc 1 l v ^f J_ iu 2 

l+4c 2 " 4c 2 4(ttJ ab ) 2 ' K ] 

thus showing that exp(tQ + .E B ) is essentially a diagonal matrix exp( y l^l^ l /T+I^I A E B ^) of 
phase factors as long as oo\ <C n\ J AB |. Noting that I^E+ still commutes with H trn , it 
follows that 

where the last two factors are conditional phases. 

It is interesting to observe that since sin 2 (| v /i+4^f') = if a/1 + 4cf = 2n for 
an integer n > 0, an exact implementation of the c-NOT is obtained when |ci| = 
|7tJ ab |/cji = \\/An 2 — 1 or uoi = 2n \ J AB \ / y/Ari 2 — 1, in which case the time required for 
the c-NOT is 7r/wi = V4n 2 - 1/|2J AB |. For n = 1, this is y/3/\2J AB \, which is longer 
than the l/\y/2J AB \ required by the on- resonance implementation, although in some 
circumstances the fact that an additional soft pulse is not needed in the on-transition 
case might be an advantage. 

5. CONCLUSIONS 

We have presented a detailed analysis of two Pound- Overhauser implementations of the 
controlled-NOT gate, one using an on-transition x-pulse with u\ <C 2n\J AB \, and the 
other using an on-resonance x-pulse of power uj\ = n\J AB \ and duration l/(y / 2|J AB |) 
followed by a soft itI A pulse. The correctness of the effective Hamiltonians for both 
implementations have been validated by NMR experiments, and the phase corrections 
derived using geometric algebra. In the course of this analysis, it was shown that the time 
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required for the on-transition implementation could be decreased to t — v / 3/(2|J AB |) 
while actually making the implementation exact (up to conditional phases) by increasing 
the power to u)\ = 2n\ J AB \/y/3. Neither the on- resonance nor the on-transition Pound- 
Overhauser implementation is as efficient as the pulse sequence in ( pory et~aT| |1998| ), 



which takes only 1/(2|J AB |) plus two soft pulses. Nevertheless, the Pound-Overhauser 
implementations have the advantage that in systems of more than two spins, it should be 
easier to find a sequence of tc pulses which refocuses the evolution of all the remaining 
spins while the gate is in progress, because it is not necessary to treat the spins to 
which the gate is being applied specially. Geometric algebra calculations with sums of 
transition Hamiltonians ( (Somaroo et al.| |1998| ) further indicate that the simultaneous 
irradiation of multiple resonances, which we call compound pulses, can be used to directly 
excite multiple quantum transitions, and thereby accomplish in a single step what might 
otherwise take considerably more time. Their effective Hamiltonians will be the subject 
of a future paper. 
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Appendix A. DIAGONALIZATION OF THE ON-RESONANCE 
HAMILTONIAN 

In the on- resonance case, the characteristic polynomial det(H c{i /uJi — A) of the effective 
Hamiltonian is a quadratic in A 2 , with eigenvalues 



A 



| (2 + c\ + c\ ± 2y/\ + <*{! + <*)) , (A.l) 



where C\ and c 2 are given by Eq. (|]). It is therefore reasonable to expect that it will be 
possible to analytically diagonalize this Hamiltonian, and we now show how this can be 
done using geometric algebra flSomaroo et al.| |1998| ) . 



We begin by rewriting the on-resonance Hamilonian (Eq. (|7])) in the form 

HeffM = 2 Cl / A / z B + I A + c' 2 J Z B exp( t/ i2/ y B ) , (A.2) 

where 



c 2 = yl + c 2 an d tan(//) = l/c 2 . (A. 3) 



Rotating spin \i about 2J B this yields 



y 

rB 



= 2ci J A exp(-t/i2/ B )/ z B + I A + c' 2 / z E 
= 2 Cl / A (J B cos(/i) + J B sin(/i)) + I A + 4J Z B 
= e^ 21 " 1 " (2 C ;/ A J Z B + 2 Cl / A / B sin( yU ) + c' 2 / B ) e 



(A.4) 
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where 



c[ = a/1 + c\ cos 2 (/i) and tan(z/) = l/(ci cos(/i)) . (A. 5) 

The outer exponentials can be eliminated by applying the opposite rotation, obtaining 



2\ + I B z e- LK + 21 * E h + - 2A_/ z B e- tK - 2/ ^ A 



where 



KiH'^/u^K = 2A +J E A I B - 2A_£ A I B 



\ IT'A jpB \ ttiA jttB \ ttiA ttiB , \ ttiA tti 

A+-G/ I Hi i — A-l-C/ I £/_ — A_-C/_.G/ i -p A_.£/_-E/ 



B 



(A.6) 



= 2^1*1* + 2 Cl / A / B sin(/i) + c' 2 / z B . 
To complete the diagonalization, this is rewritten as 

HSfM = ((ci + ^lf + dsin^)^- 

(K-c'^J^ + dsin^)/, 8 )^ (A.7) 



A± = ±vVi±4) 2 + (cisin(/i)) 2 (A.8) 
are the eigenvalues as above, and 

Ci sin(u) , , , 

tan(K± = \ V . A.9 
Ci ±c 2 

Therefore the conditional rotation 

K = e - LK+I y + e'^- 1 * E f _ 

completes the diagonalization to 



(A.10) 



(A.11; 



While it is possible to write down the time-dependent exponential of the 
diagonalized Hamiltonian, and to transform it back to that of the original Hamiltonian 
H e s, the resulting propagator is too complicated to yield much insight into the 
dynamics. Nevertheless, in the next Appendix the above diagonal form will enable us to 
give an elementary proof that off-resonance effects can be neglected in the on-resonance 
case. 



Appendix B. OFF-RESONANCE EFFECTS ON THE PROPAGATORS 

Appendix B.l. The on-resonance case 

In order to justify the above assumption that the J B term in the Hamiltonian can be 
neglected, let us extend our definition of H c s (Eq. |7|) to 

H% = 2vr J AB I A I B + (lUq - ^ B )J B + ^(J A + «/ B ) . (B.l) 
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Note that H 
by 



eff I a =0 



i3"" ff = H B g — ujiI^ as above. Accordingly, Eq. ( |A.3| ) is replaced 



\J a 2 + c\ and tan(/i) = a/c 2 ■ 



(B.2) 



The rest of the diagonalization of H^ s goes through exactly as for H cS = H eS . It now 
follows from Eq. ( |A.6| ) that exp(— LtH^ s ) 



e -i-nTf e -tvl^ g-w+E^lf e -LK-Btl* e -wiUX+E%I^ 



x e' 



(B.3) 



ye y 2 e 



where /x, i/, and A± are all functions of a. 

We wish to show that the propagators exp(— tii2"° ff ) and exp(— btH e s) are 
approximately equal. Although this is implied by first-order perturbation theory, for 
any noninfinitesimal perturbation this standard argument falls short of being a rigorous 
proof. To this end, let us fix t and define the functionjj]] 

2 



/(«) = \ 



-utH 



iff Q~ LtJrl cft 



-UH a „ aH°„ 

p err p eft 



(B.4) 



Using the fact that ||U|| = 1 for unitary U, the Cauchy-Schwarz and triangle 
inequalities, and the invariance of the scalar part under cyclic permutations, the 
magnitude of the derivative may be bounded as follows: 



\n<*)\ 



d/dae- LtH °«e UH °« 



< \\d/d> 



ae ^ff|||| e ^ff| 



\d/dae~ 



< 2\\tl°dfi/da\\ +2\\L2I~I°du/da 

->A r B o.„ /q II , oil . ET'A r B 



(B.5) 



2\\LE^dK + /da\\ +2\\iEtI*dK„/da\ 



+ \\uitiE*IZd\+/da\\ + W^tiE^dX-fda 



To bound the first two terms on the right-hand side, the derivatives are evaluated 
and simplified by means of the assumptions 2 < 1 + c\ < c|, i.e. 



2 \\il*dii/da\ 



\d/da arctan(a/c 2 ) 
1 



|c 2 |(l + (a/c 2 ) 2 ) 



< 



(B.6) 



C 2 



| Since multiplying a propagator by a scalar phase factor has no effect on how it transforms the density 
matrix, the correct way to compare two N x N propagators U\, U 2 is to compute the expression 



1 - \t r {U 1 U 2 )/N\' 



= \-{U l U 2 ) 2 -{lU 1 U 2 ) 2 

= \\\Ux- U 2 \\ 2 \\U 1 + U 2 \\ 2 + \\\Ui - lU 2 \\ 2 \\U 1 + lU 2 \\ 2 - 1 



It is sufficient, however, to show that just one of the above norms vanishes, e.g. \\Ui — U 2 \\ 2 , meaning 
that the propagators are equal without any overall phase differences, and this turns out to be the case 
for all the propagators considered in this paper. 
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2\\L2I*I*dv/da\ 



y z 



\d/da arctan(l/(ci cos(/x)))| 
8/ da arctan(A/o; 2 + c\j (cic 2 )) 
a I Ci 1 1 c 2 1 



(a 2 + c 2 (l + C 2 )) V / a T Tc 



(B.7) 



< 



Cl 



< 



1 



:i+c 2 )c 2 - 



I 12' 
I c l I c 2 



where the inequality is obtained by setting a — 1 in the numerator and a = in the 
denominator. 

To bound the terms depending on k±, define £ = a 2 + c 2 (l + c 2 ), and proceed as 
follows: 



2 ||^±Jy dre±/da| 



2 ^ 2 \d/da arctan (ci sin( / u)/(c' 1 ±c' 2 ))| 



1 

71 
i 

71 



9a 



arctan 



Cl ((c 2 



V? ± (« 2 + c 2 ) 



(B.8) 



a 



(a 2 + c 2 )((1 + a2 + c 2 + c 2 )v ^ ±20 

An elementary analysis of the denominator shows that it is nonnegative and reaches 
its minimum in the interval [0, 1] at a — in both the "+" and "— " cases, while the 
numerator is likewise nonnegative and reaches its maximum at a = in the "+" case 
and a — 1 in the "— " . This together with further simplifications lead to the bounds 



1 

71 



da 



< 



c i 



V2c 



and 



1 

71 



8k. 



da 



< 



c i 



V2\c 2 \(\c 2 \ -y/2\ Cl \) 



(B.9) 



Lastly, the terms depending on X± are 



±uj 1 tLElI°d\±/da\\ = 



UJ\t 

u>it 
471 



d/da i v / (c , i±c / 2 ) 2 + (c 1 sin(/i)) 2 



a(V?±l) 



y/£(l + a* + + ±2y/£) 



(B.10) 



As before, an elementary analysis shows that in both the "+" and "— " cases, the 
denominator is minimized at a — 0, while the numerator is maximized at a — 1. This 
together with further simplifications leads to 



UJ\t 

2^2 



dX, 



da 



< 



UJit 



and 



C 2| 



2^/2 



(9a 



< 



4(|c 2 | - y/2\ Cl \ 



Recalling that in the on-resonance case u)\ 
ui\t < ujit max = ix I y/2, we finally obtain 



ir\J I, so that ci 



1 



1 



!/'(")! < 1^ + 3 + 



+ 



+ 



7T 



+ 



tt/v/2 



l C 2| 



4 V2c 2 2 V2\c 2 \(\c 2 \ - V2) 8|c 2 | 4(|c2|-v^) 



(B.H) 
1 and 

(B.12) 
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for all < a < 1. Letting g{c2) denote the right-hand side of this inequality, it now 
follows from our assumption |c 2 | = \uq — \/u 1 ^> 1 that 



e ~ttff cff _ e ~LtH aff 



< /(0) + max 

0<a<l 

= s(c 2 ) < 1 , 



(B.13) 



as desired. 



Appendix B.2. The on-transition case 

It remains to be shown that exp(— l{tx / uji)H c q) m exp(— l{tx / uJi)H° eS ) in the on- 
transition case. This is again expected from first-order perturbation theory, but a 
rigorous proof is rendered nontrivial by the fact that H c g admits no closed-form 
diagonalization. The proof presented here is based upon the "sinch" commutator series 
expansion of the directional derivative of the matrix exponential which may be found in 
flNajfeld and Have| |1995|) , using however geometric algebra to proceed in a coordinate- 
free manner (|Somaroo et al.| |1998|) . 

Thus we define the perturbed Hamiltonian H® s = H® s + acuil^, so that 



(B.14) 



The "sinch" expansion of the direction derivative is^ 

^ {i*. aw: - ,A1 



e ~ l -2 H eS 



x ) V2 



eff 



,fe=0 



■l) k (2k 




, t rra 



where the commutator powers are defined recursively by 



{X,Y } 
{X,Y k } 



X , and 

[{X,Y k - l },Y\ (*>0). 



(B.15) 



(B.16) 



The essential thing to note is that {J x , |i3"" ff } = {I®, \H° eS }, so that the summation 
in Eq. ( |B.15[ ) is independent of a. 



% The missing factor of "t", as compared with Eq. (105) in (Najfeld and Havel 1995), is due to the 
fact that we are treating here the time multiplying the Hamiltonian and the time multiplying the 
perturbation as independent parameters, and dropping the latter. 
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In order to evaluate the commutator powers explicitly, one proceeds as follows: 

{J X B , (|/f e ° ff ) } = J X B by definition; 

{ll(lH° eS y} = c 1 |[/ x B ,(c 1 + c 2 )/ z B + 2c 1 /^ + J r x] 
= o; 1 |((c 1 + c 2 ) + 2 Cl ^)[/ x B ,/ z B ] 
= u 1 U(c l + c 2 ) + 2c 1 I*)(-±I^ • 



. 2 u^ 1 i "^i^z/v 2 y 
{II (\H%f) = ^ 1 |((c 1 + c 2 ) + 2c 1 ^)[-|J B ,| J Ff° ffJ 
m((c 1 + c 2 ) + 2 Cl /^)) 2 / x B ; 



and in general 

Ux B , (f^cVn = m((c 1+C2 )+2 Cl ^)) 2fc / x B 



It follows that the commutator series is given by 

y {g, (t^D = m((c 1 + c 2 ) + 2c 1 ^)) 2fc \ 

^(-1)*(2* + 1)! ^ (-1)*(2* + 1)! J x (B.17) 

= sinc^f (( Cl + pa) + 2c 1 I*))i x B , 

where sinc(JT) = sin(X)/X as usual. This can be rewritten in terms of scalar functions 
as follows: 



sinc^f (( Cl + c 2 ) + 2 Cl I*))(E* + Et) 

sinc(o;i|((ci + c 2 ) + c x ))E% + sinc^f ((ci + c 2 ) - c 1 ))E A . 



(B.18) 



Finally, on taking norms of both sides of Eq. ( p,15|) , one obtains 



|/ B -Ve-^|| < ||e- t l«? ff |||| s inc(a; 1 |( C2 + 2 Cl ))^I x B 



- 4 1 

+ sinc^l^^JfllUe-*^! 



< |sinc( Wl |(c 1 + 2c 2 ))|||^J x B | 
+ |sincK|c 2 )|||^ A / B || 

< ^(|c 2 + 2 Ci r + |c 2 r) 

< V2 ( 2(|c 2 | + | Cl |) \ 
" tu x \\c 2 \(\c 2 \ - 2\ Cl \) J 

< V8(i o;i(|c 2 | -2|ci|)) _1 . 

Now consider the function 



(B.19) 



/(«) = \ 



e ~ LtH eH _ e ~ UH cfi 



(B.20) 
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the magnitude of whose derivative may be bounded as follows: 



!/'(«) 



< tu^I^ ■ We~ ttH " s \\\\e UH °- s \ 

< v / 8(|c 2 |-2|c 1 |)- 1 



(B.21) 



It follows that 



/(I) 



< 



e -ttH cft _ e ~LtH { l B 



< /(0) + max \f(a) 

0<a<l 



V8 



\uj*-cu$\ -2tt|J ab | 



< 1 



(B.22) 



for all times t, as desired. 
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